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We analyze the relation between evolution equations at low x that have been derived in differ- 
ent approaches in the last several years. We show that the equation derived by Balitsky and 
Kovchegov is obtained from the Jalilian-Marian-Kovner-Leonidov-Weigert (JKLW) equation in 
the limit of small induced charge density. We argue that the higher nonlinearities resummed by 
the JKLW equation correspond, in physical terms, to the breakdown of the eikonal approximation 
when the gluon fields in the target are large. 

1 Introduction 

In recent years, there has been renewed interest in the understanding of the physics of systems with 
large number of partons. These studies have been essentially motivated by two large experimental 
programmes — low x deep inelastic scattering (DIS) at HERA and heavy ion collision experiments 
at RHIC and LHC. Both physical situations involve a large number of participating gluons. In low 
X DIS these gluons are generated in the proton light cone wave function by the evolution to low x, 
whereas in the nuclear collision this evolution is enhanced since the nuclear wave function contains 
many gluons already at moderate values of energy. 

The growth of gluon density leads to interesting physical consequences, the physical understanding 
of which has been steadily improving. One universal feature now believed to be true, is the 
saturation of gluon densities. Apparently, the number of gluons per unit phase space volume 
practically saturates and at large densities grows only very slowly (logarithmically) as a function 
of the parameter that triggers the growth. The relevant parameter could be 1/x in the low x 
regime or the atomic number of the nucleus A in heavy ion collisions. This saturation takes 



place at values of transverse momentum below a certain saturation momentum kg, which itself 
depends on 1/x and A. The nature of this dependence is less well understood. In the analysis 
based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution Q and on the double logarithmic 
approximation (DLA) ^ the dependence is powerlike kg oc while other approaches |p 

suggest a much slower dependence. In the case of power dependence, the saturation momentum 
at HERA is estimated to be in the range of 1-2 GeV with similar, slightly higher, values at LHC. 
Optimistically, one can hope that the saturation region is itself semiperturbative, that is the value 
of the coupling constant is reasonably small and, therefore, weak coupling methods can be applied 
to the quantitative analysis of the phenomenon. 

The physics of saturation must have experimental manifestations. The simplest and the most 
direct, in a way, is the unitarization of the total DIS cross section. This is, however, also the 
least interesting one. First, since the effect of unitarization is almost kinematical, one does not 
need high partonic density, it is enough to have a large number of partons, not necessarily in the 
same bin of the phase volume |5|. Second, because the experimental status of unitarization 
is unclear. So far, all DIS data on the total cross section can be reasonably well described by 
linear Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution without the need to include 
nonlinear effects @. Although physically it is hard to believe that the leading twist perturbative 
approximation can be applied at as low as 1 GeV^ and although some aspects of the gluon 
distribution that emerge form these fits are intuitively not satisfactory, present inclusive DIS 
data cannot be considered as an unambiguous confirmation of nonlinear effects. 

The realm of nonlinear effects is, however, much richer than the total cross section. In particular, 
one expects qualitative changes in the structure of the final states as one moves into the saturation 
region. The study of these effects has, however, not started in earnest yet and we have a long way 
to go before being able to make verifiable quantitative predictions. 

In particular, one needs a well defined formal framework to perform calculations. Several ap- 
proaches to the problem have been developed in recent years by different groups. They all rely 
on the smallness of the coupling constant while resumming the effects of large number of par- 
tons/partonic density. The aim of all these approaches is essentially to derive the evolution of 
the hadronic scattering cross section with 1/x. They, however, utilize different techniques and 
conceptual frameworks and the resulting evolution equations look rather different. It is the aim 
of this paper to explore the relation between some of these different approaches in an attempt to 
understand where they diverge from each other in terms of physics input. 

In particular, we will concern ourselves with three recent works Ref. 0, Ref. and Refs. ||9-13|. 
In Ref. the evolution equation for the scattering amplitude is derived using the effective action 
and the eikonal approximation in the target rest frame. Ref. uses the dipole model method 
of Refs. |14, 15|. And, finally, Refs. P-]l^ uses the effective action in the projectile rest frame to 



derive the evolution of the hadron light cone wave function with 1/x. We will refer to the resulting 
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evolution equation as the Jalilian-Marian-Kovner-Leonidov-Weigert (JKLW) equation. 

The outhne of this paper is the following. In Sec. |2| we rederive the evolution equation of Ref. in 
a simple and intuitive way. This derivation makes it obvious that this approach is equivalent to the 
approach of Ref. up to subleading corrections in l/Nc- This is not new and was noted already 
in Ref. ||8|. In the following, we will refer to this evolution equation as the Balitsky-Kovchegov 
(BK) equation^. We discuss the physical picture of this evolution and resulting unitarization of 
the total cross section in both target and projectile rest frames and point out the effects due to 
which the approximations involved should break down at extremely small x. The breakdown of 
the approximation should have very little effect on the unitarization of the total cross section, 
since, especially for large targets like nuclei, the black disk limit should be reached while the 
approximation is still valid. However, one does expect the structure of the final states to be 
strongly affected. Our discussion here is, in large measure, parallel to that of Refs. 

In Sec. ^ we relate explicitly the calculation of Ref. to that of Refs. In particular, we 

calculate the basic physical quantities appearing in the evolution equation of Refs. p|-[T^ in the 
approximation of Ref. 0]. We show that the results of Ref. Q are recovered from Refs. ||9-12] in 
the limit of small induced fields. We also show that the double logarithmic limit of the evolution 
of Ref. is trivial. That is, in the double logarithmic limit, the evolution equation for the gluon 
distribution function (defined operatorially as the number of gluons in the lightcone gauge in the 
infinite momentum frame) becomes linear and does not contain any Gribov-Levin-Ryskin (GLR) 
type corrections^. This is in contrast with the result of Ref. where it was shown that the 
double logarithmic limit of the evolution Refs. [p[-[T^ results in a nonlinear equation. We point out 
that this is indeed a very natural result from the point of view of the dipole model approach. 

In Sec. § we transform the full calculation of Refs. [^12| into the framework of Ref. 0. We 
show that in the approach of Ref. it corresponds to abandoning the eikonal approximation or, 
equivalently, to the inability to fully describe the target by a classical field. We also point out 
technical reasons which lead us to believe that, in fact, in the framework of the effective action of 
Ref. such a failure is expected when the evolution is continued to very low values of x. 

Finally, we conclude with a brief discussion in Sec. ^. 



^The 0{1/Nc) difFerences between the equations derived in Ref. and Q do not carry essential new physics. 
They therefore do not affect our understanding of the relationship between the generic frameworks of the BK and 
JKLW equations. 

^This is not to say that the evolution of the DIS cross section for which the equation of Ref. Q and ^ has been 
derived is linear in the double logarithmic limit. The GLR type nonlinearity does indeed appear in the evolution 
equation for the virtual photon cross section due to the nonlinear relation between the cross section and the gluon 
distribution function. 
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2 A simple derivation of the BK equation 



In this section we will give a simple derivation of the evolution equation first derived in Ref. |^ 
and discuss the physical picture behind it. Consider the deep inelastic scattering at low x. We 
will work in the frame in which the photon fluctuates into an energetic quark-antiquark pair long 
before it reaches the target, but where most of the energy resides in the target hadron which moves 
very fast. The scattering of the quark-antiquark pair is dominated by its interaction with the 
gluons in the target. Since the target hadron moves fast, the time evolution of the gluon fields 
is slowed by Lorentz time dilation. Also, due to Lorentz contraction, the gluon fields are well 
localized in the plane perpendicular to the direction of motion, which we take to be the positive 
X3 axis. The target can, therefore, be modeled by a distribution of static gluon fields localized at 
x~ = 0. As the scattering energy increases (x decreases) the gluon fields of the target change due 
to contributions of quantum fluctuations. It is this evolution in x of the hadronic ensemble that 
we intend to describe in terms of the evolution equation. 

2.1 The BK equation. 

In this section we will use the lightcone gauge A~ = 0. In this gauge, following Ref. Q, we take 
the vector potentials representing the relevant gluon field configurations to be of the form 

= 0, 6+ = b{x±)6{x') . (1) 

Here and in the rest of this section, unless otherwise specified, we use the matrix notation for the 
gauge field b'^ = b^ta etc., where ta are the generators of the SU{N) group in the fundamental 
representation. One can reasonably ask whether the vector potential of this form is the only 
relevant one. This turns out to be a nontrivial question. In fact, we will argue later in Sec. I that 
this is not quite the case if we want to be able to describe the evolution up to arbitrarily small 
values of x. At this point, however, we follow Ref. [^. We will return to this question in Sec. ^ 

The DIS structure function can be written in the following general form 

F,ix,Q') = -^ [ ^^^^^^x^-y^,z)N{x^,y^,y). (2) 

Here, x± and y± are the transverse coordinates of the quark and the antiquark in the pair, z is 
the fraction of the pairs longitudinal momentum carried by the quark and y is the rapidity of 
the slowest particle in the pair. Also, ^{x± — y±,z) is the square of the "wave function" of the 
photon — the probability that the virtual photon fiuctuates into the pair with given coordinates 
and momenta — and N{x±,y±,y) is the cross section for the scattering of the pair. 

The wave function $ is well known. It is given for example in Ref. |8|, but its explicit form will 
not be of interest to us. We concentrate our discussion on the scattering cross section N. If the 
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quark-antiquark pair is energetic enough the scattering cross section is eikonal 



N{x^,y^)=tv{V{x^)V^{y^)-l)^, (3) 
where V (V^) is the eikonal phase for the scattering of the energetic quark (antiquark) 



V{x~^ = 0,x±) = V ex.p — ig dx A'^ {x'^ = 0,x±,x ) 



(4) 



with the vector potential in the fundamental representation. We, therefore, have to calculate the 
average of (y{x~^ = 0,x±)V'^ {y^ = 0,y±))^ over the hadronic wave function as indicated by 
(. . . )a. In our frame, the quark and the antiquark move with velocity of light in the negative X3 
direction. All the fields in Eq. (^), therefore, have vanishing x^ coordinate. This will also be the 
case for all the fields in the rest of this section. For simplicity, we suppress the x~^ coordinate in 
the following. 

In the leading approximation the vector potential is given by Eq. (1^) and the scattering amplitude 



IS 



N{x^,y^) = (tr[U{x^)UH 



y±) 



with 



+00 

U{x±) = 'Pex.p —i J dx~b'^{x±,x~) 



(5) 



(6) 



To calculate the order as correction to this expression we write the vector potential as 

A+ = -b+ + a+ 
9 

with being a small fluctuation and expand the eikonal factors to second order in 



(7) 
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Recalling that the classical background vector potential is a delta function in x , we have 

+00 

ig I dx a' [x±,x 



-i-00 

V{x±)=Vexp — ig J dx~a~^{x±,x^) U{x±)'Pexp — ig Jdx~a'^{3 

-00 
+00 

= Uix,)-^g{ ldx-a-ix,,x~Mx,) + Uix,) J dx-aHx,,x^)} 

-00 


~9^{ j dx^dy'9{y' - x~)a'^{x±,x~)a~^{x±,y~)U{x±_) 
—00 

+00 
+ dx^a'^{x±,x~)U{x±) dx~a^{x±,x~) 



+00 



+ Uix,) Jdx-dy-Oiy- - .-)a+(x„ .-)a+(x„ ,-)} . 


All contributions break down into x~ ordered pieces because of the x~ structure in Eq. ([l|). 
Now, together with the analogous expansion for VK we insert this into Eq. (^) and obtain 

tr{Vix^)V\y^))^ - tr{U{x^)U\y^))^ = 

= £/^tr// dw'a+{xj_,w')U{xj_)U^y±) dz'a+{yj_,z~) 
J— 00 J —00 

/O r+00 
dw~a'^{x_i,w^)U{x±) / dz^a'^{y±,z^)U\y±) 
-00 Jo 

/■+00 t-O 

+ U{x±) / dw~ a'^ {x±,w^)U^yj_) / dz~a~^{yj_,z~) 
Jo J -00 

r+oo r+00 

+ U{x±) / dw^ a'^ {x±,w^) / dz^a'^{y±,z^)U'^{y±) 
Jo Jo 



-00 




dw dz 9{z — w )a'^{x±,w )a'^{xj_,z )U{x±)U'^ {y±) 

o 

r+00 

dw~ a'^ {xj_,w^)U{xj_) / dz~ {x±, z^)U'' {y±) 

o Jo 
r+00 

— U{x_l) / d'w~dz~6{z~ — w^)a'^{x±,w~)a'^{x_L,z~)U\y_L) 

Jo 

— U {x j_)U'^ {y j_) / dw^ dz^ 9{w~ — z^)a'^ {y±,w~)a~^ {y±, z~) 

J —00 

r+oo rO 

-U{x±) dw'd^{y±,w')U\yA_) I dz' a'^ {y±, z') 
Jo J -00 

t-+oo 

-U{x±) dw~dz~e{w' - z')a~^{y±,w')a'^{y±,z')U'^{y±)) . 

In I b,a 



In writing Eq. (^) we have anticipated that {a^)^ = as in the free case. This can easily be shown 
using the exphcit expression for the fluctuation propagator given below. 



Although this expression is a little cumbersome, the physical meaning of the various terms is very 
clear. Put more compactly the structure of the above is determined by 



y±) 
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Diagrammatically the r.h.s. can be represented as follows 



6bZ ' ^'5bt 
-U\yi_))) 



(10) 




(11) 



The straight lines represent the eikonal factors f7, while the curly lines denote the gluon fluctuation 
propagator (a^a+)^, evaluated in the fixed background The terms in Eq. ([lO|) with first 

order derivatives correspond to processes where the gluon is emitted by the quark and absorbed 
by the antiquark (or vice versa). Those will be hereafter referred to as "exchange contributions". 

The terms in Eq. (|lO|) with second derivatives acting on U (or C/^) correspond to the diagramms 
where where the quark (or antiquark) emits the gluon and then reabsorbs it at a later time; typical 
self energy corrections. To contrast them against the exchange contributions we will also refer to 
them as "non-exchange contributions". 



In Eq. (^) we are looking at a x ordered breakdown of the diagrams in Eq. (jl^) with the first four 
terms summing up to the exchange contribution and the rest to the non-exchange contributions. 

With the vertices known, we only lack an explicit expression for (|a+a+)^. The QCD action 
expanded to second order in the fluctuation field a"^, in the presence of the classical background 
6"^, in our lightcone gauge is 



S = \[<[- {d'fK - 2(5^a+)(a-<) 



(12) 



Recall that we are interested only in the propagator of the fields at equal x^. Consequently, it is 
only the on-shell part of the propagator that is relevant for our purposes. We can, therefore, use 
the classical equation of motion for a"*" 



d- 



(13) 
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Substituting this in Eq. ( [T^ ) we get 



S = --aUD\b6'^ai, 



where 



2D+[b]d--{d^n 



•ab 



(14) 



(15) 



The propagator of the spatial components of the vector potential is, therefore, simply given by 
—i/D'^. The explicit form is very simple and can be found, for example, in Ref. |16|. 



1 






-j 

ab J 



dp 



2p-(27r)3 



[e[x- - y-)9{p-) - e{y- - x-)e{-p-)] X 



(fp±(fq^e''P-''+''^-y 



(27r): 



(16) 



with p~^ = ^p, = and q = p ■ The adjoint color matrix^ 

U-,Hx~,y-,z±) = {eix-)e{y-) + 9{-x-)9i-y~))5ab 



+ e{-x-)9{y-)Uab{z^) + e{x-)e{-y~)ul,{z^) . 



(17) 



represents a phase factor one picks up when crossing the x~ = plane due to interaction with a 
field of type Eq. (|T|). Here, Uab{z±) is the adjoint version of the fundamental U in Eq. (^). If the 
x~ = plane is not crossed the propagation remains free. 

We can now write the on-shell correlator of the "+" component of the vector potential as 

(a+(x+ = 0,x±,x')a^{y^ = 0,y±,y~))^ 

9* • _ • _ 

= {7^aa{x+ = 0,x±,x )al(y+ = 0,y±,y )-^) 



-didi I ^^—[e{x- - y-)e{p-) - e{y- - x-)e{-p-)] 



(p-)3 47r 



'e 2p- 



2p- 



(2^)2 (2^)2 

X {e{x-)e{y-) + e{-x-)9{-y-) 

+ e{-x-)e{y-)U{z^) + e{x-)9{-y-)W{z^)} . (18) 

This expression displays a separation into x~ ordered contributions that seamlessly matches up 
with what we have already seen for the vertices in Eq. (P). Diagrammatically, the x~ ordered 



■^More rigorously, the structure of U^^ is given by U~i^{z±) = e*'*'^ )s(.y ))>'(^±) ^ However, the difference between 
this expression and that given in Eq.(p^ only shows up if it is multiplied by derivatives or S(x~) factors. Since 
we encounter no such factors in our calculation we will be using Eq.(jl^) throughout. 
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exchange contributions are given by 





while the (connected parts of the) non-exchange ones are represented by 



(19) 




(20a) 



(20b) 



Whenever the x =0 plane cuts the fluctuation propagator, the diagram contains a factor , 
otherwise the fluctuation propagator is free. 

Algebraically, the corrections to the scattering cross section involve integrals of this propagator 
with respect to x~ and from zero to either + or — inflnity. This is straightforward to do. For 
example, for the second term in Eq. (ph we need 



dw J dz {a'^{x±,w )a'^{y±,z )) = 



+ 00 

dn^ I f o f d'^m d'^a\ p*Px(a;±-2±)+jg±(z±-s/±) _ 

dldl I / d'z^ I U^biz^) 







p vr 



(2^)2 (2^)2 



+ 00 

1 fdp~ d\~ d\ 

IT J p- dl dj_ 





(21) 



Here, {x±\0\y±) means the matrix element of the operator O in the coordinate basis in the usual 
sense. We treat U as an operator in the coordinate space with matrix elements {x±\U\y±) = 
U{x±)5{x± — y±) and the products in the last line are understood in the operatorial sense. Ex- 
plicitly, 



+00 



Jdw j dz (a+(x_L,'u; )a^{yj^,z )) = ^ ^ 



•i-K^ J p 




d^z ^-'^t;^~^^ -U^(z^- (22) 



(x-z)i(y-z)i 
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In the same way we obtain for the other contributions with interaction with the background 

+00 

-00 

T^J Q P dj_ ""dl 47r'^J P J {x - z)j_{y - zyj_ 

+00 

- J dw~ J dz~{a^{x±,w~)a^{x±,z~)) = 
-00 

= -- r^^{x±\tUa,^\x±) = I d'z^^^ ■ Ua,{z^) (24) 

TTj P di di 47r'^J P J {x-z)j_ 



+00 

- J dw~ J dz'{a^{y^,w~)a^{y±,z~)) = 
-00 

To simpUfy the color structure of these expressions we use the identity 

Uab{z±){taU{x^)r''{tbUHys.)y' = 2tv[taU{z^)tbU\zi^)] {taU{xs.)r''itbU^iy±)r' 

\n4u{z^) ■ uHy±)r\u\z^) ■ uix^)y^ - u{x^rfuHy±r'] . (26) 



2N, 



Note that the integral over the frequency logarithmically diverges. In fact, we have to integrate 
only over a finite interval of frequencies. The gluon field modes of very low frequency have been 
already included in the background field 6+ and, therefore, the fluctuation fields at these low 
frequencies should not be considered. The lower cutoff on the frequency of the modes that are 
being integrated is inversely proportional to the initial value of xq at which we start the evolution. 
The upper limit on the high frequency side is furnished by the maximal rapidity of the quark (or 
antiquark) in the virtual photon which is of the order of 1/x. The ratio of these two cutoffs is of 
order xq/x. Thus, in the leading logarithmic approximation we identify 



/ 



^ = ln^. (27) 

p~ X 



The calculation of the remaining contributions (the ones with no interaction with the background) 
proceeds along similar lines and is given in the appendix. 
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Collecting all contributions together we obtain 



— Inf— ~ 



tr(y(xx)Ft(yx)>^ - tr(C/(xx)t/t(yx)> 

2tr([/(xi)C/t(z_L))tr([/t(yj„)C/(z^)) - 2iV,tr(C/(x^)C/t(y_|_)) 

ti{U^{y±)U(z^))ti{U\z^)U{x^)) - NMuHy±)U(x^)) 



{x - z)±-{y - z)± 

{x-z)l{y-z)l 
1 



(x - z)l 

1 



{y-zY 



(28) 



The eikonal factors V themselves should be considered as functions of x, so that U = V{xq). Differ- 
entiating this equation with respect to In 1/x we recover the evolution step for tr{V{x±)V^y±))A- 
This is precisely what was found in Ref. 0^ 



At large Nc the products of traces in Eq. ( pS] ) factorize: 

(tr([/(xx) . uHz±))tT{uHy^) ■ U{z^))) {ti{U{x^) ■ uHz±))){ti{U\y^) ■ U{z^)). (29) 

Equation ( |28| ) then becomes a closed equation for the evolution of A^(x_l, yj,) = (ti{U{x)U\y) — 1)). 
It is identical to the nonlinear evolution equation of Ref. p] . 



2.2 The physical interpretation. 

Now, let us discuss the physical picture of this evolution. As always with DIS, the physical picture 
depends on the frame in which one chooses to view the process. We have specified the frame to 
some extent by declaring that the photon fluctuates into a qq pair long before the target. However, 
we are still free to put the subsequent evolution in x either into the evolution of the photon wave 
function or into the evolution of the gluon field distribution in the target. We will refer to the former 
picture as the "projectile evolution picture" and to the latter as the "target evolution picture". 
The calculation is, fortunately, noncommittal on this point and we will consider both pictures in 
turn. 

In the projectile evolution picture, the higher energy of the scattering is achieved by boosting the 
qq pair. In this picture, the quark and antiquark have very high energy and consequently their 
wave function develops extra gluon components. The growth of the cross section with 1/x then 
is interpreted as due to the scattering of extra gluons in the projectiles wavefunction. This is 



precisely how the low x evolution is viewed in the dipole model of Mueller |14, 15|. Our calculation 



of this section has a simple interpretation from this point of view. The quark and the antiquark 



As it stands, Eq. ( |28[ ) does not provide a closed equation — it has to be supplemented by evolution equations 
for arbitrary products (V/^' ® . . . (8) kS^')a- The evolution of these higher correlators are derived following the same 
procedure as described above and leads to the full set of operator equations derived in Ref. We will give a 
compact representation of the whole set of the evolution equations in Eq. (^^, Sec. ^ 
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is the pair of pointlike color charges moving with velocity of light and located at = and 
transverse coordinates x± and y±. These color charges carry with them a fluctuating gluon field. 
When the pair is boosted to higher rapidity the gluon fields "freeze" due to time the dilation 
and become static. In the approximation when the gluon fields are frozen, they are given by the 
Weiszacker- Williams static (p^ = 0) fields created by the qq pair. In the leading order in the 
Weiszacker- Williams (WW) fields are small and are emitted independently by the quark and the 
antiquark. The total WW field is 

A' = 9^^[j:^+k]^ (30) 

where j~ {jq) is the color current due to the quark (antiquark) which in our frame has only a 
"— " component. For pointlike quark and antiquark the charge densities are delta functions in the 
transverse coordinates and in The WW field is therefore 



A\p-,z^) = g—[T- ^ + t'^ ^] . (31) 

p [x — zy [y — zY 



Here, r and r' are fundamental color matrices corresponding to the orientation of the quark and 
antiquark wavefunctions in the color space. Their exact form does not matter for our purposes. 
The WW field, if written in the particle basis, can be thought of as representing equivalent gluons. 
The number of gluons at a given transverse position is given by the familiar expression 

/•oo 

nww{zi_) (X I dp'p' trF^'ip' ,z±)F~\p' ,z±) 

°° dp- 3f-z^ jf-z^ , y'-z\ 

p- ^'^^{x-zY^^ {y-zf^^{x-zr^^ {y-zf- ^"^ ^ 







If we do not take the trace over the color indices, this expression gives the probability to have 
one extra WW gluon in the wave function of the qq pair (at the transverse position z± with a 
particular color orientation). These WW gluons scatter on the gluon field of the target eikonally 
just like the quark and the antiquark, apart from the fact that they carry adjoint charge, and so 
their eikonal amplitude is given by U rather than U. The terms in this expression are in one-to-one 
correspondence with the real contributions in Eq. (^ , that is the terms in which the gluons interact 
with the target (background). The rest of the terms in Eq.(^ - the virtual terms - as usual serve 
to restore the correct normalization of the wave function. 

To summarize, in the projectile evolution picture our calculation describes emission of the WW 
gluons into the wave function of qq long before the scattering. The transverse coordinates of these 
gluons are frozen due to the Lorentz time dilation. Subsequently, both the q and q, and also the 
gluons, scatter eikonally and independently of each other on the target gluon field. Clearly, this 
picture is identical to the dipole evolution picture of Mueller which was used in Ref. [^] to derive a 
nonlinear evolution equation. The only difference is that the dipole model uses the simplifications 
in the color algebra which arise in the large Nc limit. 
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The calculation presented above also has a simple interpretation in the target evolution picture. 
In this picture, it is the target rather than the projectile that is boosted when going to lower 
X. As already made explicit above by writing (. . . )b, one should think about the target as being 
represented by an ensemble of the configurations of . The corresponding statistical weight Z[b] is 
determined, of course, by the structure of the target at the relevant resolution scale. We will have 
more to say about it in Sec. ^. The boost of the target freezes the gluon field fluctuations around the 
target background 6"*" and, consequently, some field modes which were not important at higher x are 
now capable of inducing scattering. Thus, the ensemble of the relevant field configurations which 
characterizes the target changes. In fact, every 6"*^ now forks into a "subensemble" 6^' = 6^ + a^. 
In the weak coupling regime o"*" have Gaussian distribution with the width determined by the 
inverse of their correlation function Eq. ([I^). One can work back from here and calculate the 
modification of the distribution of the background fields. We will do this in the following sections. 

The fiuctuations of a"*" are, therefore, considered in the target evolution picture as modifying the 
ensemble of the target background fields very much like in the approach of Refs. [P-|l^,|T7|. 



2.3 Unitarization in different approximations. 



From what has been said so far, it is clear that although the calculation presented in this section 
includes into the evolution some nonlinear effects, it is not the end of the story. At very low x, 
this approximation should break down. There are clear reasons why this should happen in both 
pictures. In the projectile evolution picture, it is not true indefinitely that the WW fields are 
emitted independently from the partons in the projectile. Due to the evolution, more and more 
gluons are emitted into the wave function of the projectile and so the density of partons grows. 
At some point, the approximation of independent emissions as well as of independent scattering of 
the partons on the target must break down. This is the point at which, in the parlance of Ref. 
the pomeron loop diagrams must come into play. In the target evolution picture, the problematic 
point is the eikonal approximation for the scattering of the qq pair. At not very low x the target 
fields are not too strong. Since the qq pair is very energetic, the eikonal approximation is perfectly 
valid. However, with the evolution the strength of the target fields grows. The energy of the qq 
pair, on the other hand, stays fixed. When the fields are strong enough the quark and antiquark 
will start losing a finite fraction of their energy and, therefore, the no recoil eikonal approximation 
cannot stay valid indefinitely. 



The effect of the nonlinear evolution Eq. ( |28D on the behavior of the total cross section was studied 
in Refs. IS]. It was concluded that the nonlinearities slow down the BFKL type rise of the 
cross section and lead to its unitarization so that the cross section approaches the black disk 
limit. We want to conclude this section with a comment on the nature of the unitarization in 
this approximation. Essentially, the unitarization is brought about by purely kinematical effects. 
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This is especially clear in the projectile evolution picture. At the initial value of x = xq one starts 
with the qq pair as the only relevant component of the photon wave function, which has a certain 
probability Pqq to scatter on the target. So, initially, the total scattering probability Px^ is 

Pxo = Pqq ■ (33) 

At lower x = xq — 5x the wave function also contains a component with an extra gluon. Let the 
probability to have an extra gluon in the wave function be A and the probability for this gluon 
to scatter on the target Pg. In the linear approximation (the BFKL limit) the total probability of 
scattering is additive 

= (1 - A)Pg5 + ^{Pqq + Pg) = Pgg + A Pg . (34) 

However, this is, in fact, overcounting, since there are events where both the gluon and the qq pair 
undergo scattering and those events are counted twice in the linear approximation. One should, 
therefore, subtract the probability of these double scattering events from the total probability. This 
deficiency is corrected by writing 

P, = (1 - A)Pgg + A{Pgg +Pg- PggPg) = Pgg + A{Pg - PggPg) = + + A(l - Px + 5x)Pg ■ 

(35) 

At arbitrary low x, the same argument leads to a similar expression where Px+Sx denotes the total 
scattering probability of the projectile (which itself contains the qq pair and some number of gluons) 
at a slightly higher value of x. This is precisely the nonlinear term in the evolution equation Eq. (|28[) 
with the only difference that the extra gluon in the wave function can have arbitrary transverse 
coordinate and one should, of course, integrate over this extra degree of freedom. It is clear that 
this negative nonlinear correction leads to the unitarization of the cross section since as Px tends to 
unity the emission of the extra gluon does not increase the total scattering probability. This effect is 
somewhat similar to the Glauber mechanism, not in the sense that each parton undergoes multiple 
scattering, but that the unitarization is of a purely geometrical nature. A similar discussion, in 
the framework of the dipole model, is given in Refs. 0,^. 

In the next section we will show how to relate the approach just discussed with that of Refs. |9|- 
l|,|l|. 



3 The JKLW equation and the small induced field limit. 

We start this section by recalling the framework and results of Refs. (9-13, 1^. 
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3.1 The JKLW equation. 

In this approach, following Refs. ||20|, the averages of gluonic observables in a hadron are 
calculated via the following path integral^ 

{0{A)) = j Da'DM'0{A)Z[a] x 

X exp I - i j d^x^tiFf'Tf,^ ~W j ' (36) 

where the gluon field strength tensor is given by 

Fr = O'^A: - - gUcA'^A'^, (37) 
and W is the Wilson line in the adjoint representation along the x"*" axis 

+ ^9 j dx^A~{x"^,x^,x±)Ta . (38) 



The hadron is represented by an ensemble of chromoelectric fields, localized in the plane x =0, 
of the form 

f^' = ^Six~)a\x^), (39) 

where the two dimensional vector potential a^{x±) is "pure gauge" 

d'ai - - fahAoii = . (40) 
In Eq. (^), Z[a\ is the statistical weight of a configuration ai{x±) in the hadronic ensemble. 



The evolution in Refs. p|-p]^| is derived in the target evolution picture where decreasing x corre- 
sponds to boosting the hadronic target. This leads to freezing of part of the gluonic degrees of 
freedom. Integrating out these slow modes of the vector potential generates the renormalization 
group equation, which has the form of the evolution equation for the statistical weight Z ||9-12, 17| 



d rl S 

-Z = as{- 



2 



dln^ \2 6a{u)6a{v) ' 6a{u 

In the compact notation used in Eq. (|4l|), both u and v stand for color and rotational index and 
transverse coordinates, with summation and integration over repeated occurrences implied. This 
evolution equation for the statistical weight can be rewritten as the set of the evolution equations 



^An alternative form of the effective action was suggested in |2^] where it was also shown that it leads to the 
same evolution equation. 
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for the correlation functions of the chromoelectric field 



din 



-{ai\ixi)...ailixn)) 



a. 



Yl < < (^i)-«^'r-i {xi-i)a]i+^\ ixi+i)...ail {xn)a% (xi)) 

.0<«<n+l 

+ Yl < ai\(a;i)---«ir-i(^--i)a^™+u^™+i)-" 

0<in<k<n+l 



X aa\.\ {Xk^lWalW {Xk+l)...all {Xn)XaZX (^™' ^fc)) 



• (42) 



The quantities x[q^] and a[a] have the meaning of the mean fluctuation and the average value of 
the induced vector potential which arises from the field modes which become frozen due to extra 
boost of the hadronic target. In the leading logarithmic approximation of Refs. [^12| the two 
quantities x s^-nd a completely specify the low x evolution. We give here an explicit expression for 
the mean fluctuation x which will be the focus of our interest throughout this section 



2{x^\{^[D 



s- 



\ — 



}ah\yi.) ■ 



(43) 



For convenience, we have defined 



a 



ab 



fabca. 



I 

c ' 



Dab = 9' Sab + Olai, . 



(44) 



The operator S in Eq. (43) is given by 



— ] 



(45) 



3.2 Where does it come from? 



Technically, these results are derived as follows. One considers the quantum corrections in the 
classical background field Eq. (^). The calculation is performed in the lightcone gauge A'^ = 
with the residual gauge fixing d^A'^{x~ — > — oo) = which fixes the gauge completely. In this gauge 
the chromoelectric field Eq. ( |39| ) corresponds to the background vector potential 

b' = 9{x-)a'{x±_). (46) 

Note that, as opposed to the previous section, here we are using a different lightcone gauge: A'^ = 0. 
As a consequence, the background vector potential has a different form. 
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The complete set of on-shell small fluctuation solutions of the classical equations is 



%-,r = e 



ip 



-X } exp I I- — X 



2p- 



+ 9{x )U{x±)exp{i X 

2p 



ip±x± 



{p±) + e{x-)j\ 



(47) 



Here, r is the degeneracy label, which labels independent solutions with the frequency p^ . In the 
free case it is conventionally chosen as the transverse momentum, {r} = {p±}. The matrix U{x±) 
is the SU{N) matrix that parameterizes the two dimensional "pure gauge" vector potential a*(xj_) 



a {x±) 



iU{x±)d'uH 



The auxiliary functions 7+,^^ are all determined in terms of one vector function. Choosing this 
independent function as vl_ we have 



2D 



fjk _ pk 

tjk _ pk 



(48) 
(49) 



where we have defined the projection operators 

= i5*J 



Dl ' 



d 



,2 ' 



QiQi 



dl 



(50) 



The 7+ piece of the eigenfunction Eq. ( [4 71) is responsible for the induced vector potential since this 
is the only contribution that does not vanish at x~ — > oo, so that 

(51) 



Xafe(a;±,y±) = 47r / dp (7+,„(a;_L,p )7+,b(y±,-P )) • 



Note that the essential nonlinearity of the expression Eq. (43) is due to the denominator in the 
operator S"^^ Eq. (|45|). The reason this arises is due to the nontrivial normalization of the small 
fluctuation eigenfunctions. As discussed in detail in Refs. |^12| the proper normalization of the 
eigenfunctions requires u !_ to be chosen as a complete set of eigenfunctions of the two dimensional 
Hermitian operator 



|((_i)0-'((-i)fiO''±.!/±) = {^xKi 



D 



\y±) , 



ab 



such that 



/ 



d r^v'_ {xi_)v^ Ayi_) 



AttIp 



(52) 



(53) 
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This nontrivial normalization is the consequence of the presence of the 7+ piece in the solution 
Eq. (|47|). Equation (^T]), supplemented by Eq. (^) and the normalization Eq. (53), leads to the 
final expression Eq. (HSl). 



If the contribution of 7+ could be neglected in the normalization condition, the normalization of 
the eigenfunctions would be trivial and we would have O = 1 in Eq. (^3|). One can consider the 
limit in which 7+, or equivalently X) is small. In the leading order in the expansion in 7+ we have 
a very simple expression for x 

yS^.,y.} = 4(.xl[«'{^ + \ - ifyD.^^ - ^D,8,l,)D']jy,) . (54) 

Note that this is a different limit than the one in which the JKLW evolution reduces to the BFKL 
equation [pf-p!^. The BFKL limit corresponds to the expansion in powers of the background field 
a*. Now, we are not assuming that a* is small, but rather that the correction induced by the 
evolution, 7* is small. 



3.3 BK to JKLW: transforming between the gauges. 

We will now see that Eq. (U) is reproduced precisely by translating the calculation of the previous 
section into the language of the JKLW evolution. 

In the previous section, following Ref. 0, we used the gauge A~ = 0. This is a very convenient 
gauge from the point of view of the projectile evolution since the eikonal amplitudes in this gauge 
are given by simple Wilson line factors. We will refer to this gauge as the "projectile lightcone 
gauge", or the "projectile gauge" for short. The JKLW approach on the other hand uses the 
= gauge, which is convenient for the target evolution picture since it simplifies the relation 
between the distribution functions and the correlators of the gluon fields. We will call this gauge 
the "target lightcone gauge" or, simply, the "target gauge". Our immediate aim is, therefore, to 
calculate x using the results of the calculation in the projectile gauge. 

To do this note that the relation between the fields in the target and projectile gauges is given by 

-B^' + A^ = Vi-b^" + a^)W + -Vd^^V^ . (55) 
9 9 9 

To simplify the notation, from now on we will denote the fields in the target {A'^ = 0) gauge by 

capital letters and fields in the projectile (a~ = 0) gauge by lower case letters. This we do for both 

the background part of the field and for the small fluctuation part. The field dependent matrix V 

is given by 



V = Pexp 



dx (6+ + 50+^ 

00 



(56) 



®If so desired this expression can be written in a simple form in terms of the unitary matrix [/, since operatoriaUy 
— Ud^W . In Fourier space this gives convolutions of U{p) and powers of transverse momentum. 
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The condition = does not by itself specify the lower limit of the integration over x~ in the 
exponential. However, choosing this limit to be at minus infinity ensures that V{x~ — oo) = 1 
and, as a consequence, A'^(x~ — > —oo) = a^{x~ — oo). The projectile gauge fields satisfy the 
standard vanishing boundary conditions at infinity. This choice of the lower limit of the integration, 
therefore, guarantees that the target gauge fields also vanish at x~ — > — oo and, further, satisfy the 
residual gauge condition d^A^{x~ — > — oo) = that was imposed in Refs. |9-12|. To calculate x we 
only need to consider the linearized relation between the small fluctuations of the fields in the two 
gauges. To do this we need to expand V to first order in a^. This has been done in the previous 
section. Taking only linear terms in in Eq. ( ^a| ) and substituting them into Eq. (^) we find for 
the transverse components of the field 



4(x) = e{-x- 



f 



dx d^a 



(•0 i-x 

+ 9{x-) C7"^4(a;) - ^a6( / dx-a-^ + Ubc / dx-a+) 

D Jo 



ab 



( W)fe . 



(57) 



Here, the matrix U is the same as in the previous section and is related to the classical background 
by 



(58) 



U{x^) = V / dx b }, 

= e{x')iUd'W 



as per Eq. (|^). We have also defined 
W{x~) = Pexp 



dx-b-^ix-) } = e{-x-)l + 9{x-)U , 



(59) 



which is essentially the classical part of V . The operator ^ in the last line of Eq. (|57| ) is defined as 
the integral from — oo[[[ We will further simplify this expression by using the on-shellness conditions 



{2d-D+[b+]-dl)a' = Q. 
The resulting relation between the on-shell transverse fields in the two gauges is 



(60) 



A'^ = e{-x-){t-ira?,{x-) + 

+ e{x-)[lJab{t-ira^,{x-) 



2D, 



ab 



L92 



al{x' 



0~) - C/fe, 



92 



0^ 



}. (61) 



^ Eq. has been derived also in Refs. |^,^. The only difference in our derivation is that the meaning of the 
pole is entirely unambiguous and, as discussed above, is dictated by the residual gauge condition. 
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Note that we have to specify on which side of = the fields are taken since the solutions of 
the small fluctuation equations in the projectile gauge are discontinuous at zero. Now, recall that 
a* satisfies, everywhere except at x~ = 0, the free equations of motion. With this in mind we can 



compare this equation with Eq. (47). We see that Eq. (61) is indeed precisely of the form Eq. (47) 
with 



Y+,a = -2^afe ^ 0-) - Ubc-^aiix- ^ 0^ 



Uab{t-lfai{x~ ^0+). (62) 



Remembering that (see for example Ref. |16| ) 

a\x- ^0+) = Wa\x- ^0-), (63) 
we see that the functions v!^ , vl_ and 7!^ are related precisely by the same relations as in Eq. 



We have established, therefore, that if a* satisfies the equation of motion in the projectile gauge, 
then the transformed field of Eq. (55) satisfies the equations of motion in the target gauge. 



The only remaining question is that of the normalization of the eigenfunctions. Recall that the 
functions a* in the calculation of Ref. which was reproduced in the previous section, were 
normalized in the same way as the eigenfunctions of the free theory. That is to say, the full set 
of on-shell eigenfunctions is obtained by choosing ap^{x~ ^ 0~) as a complete set of normalized 
eigenfunctions of the unit operator in the transverse space 

/ ^ 0-)4*Jy±,x- ^ 0-) = S'^6{x^ - yx) • (64) 

Since t — / is a unitary operator, Eq. ( |62| ) tells us that vl_ is also normalized to unity rather than 
to a nontrivial operator O as in the JKLW calculation Eq. ([5^). Using this, as well as the relations 
Eq. ( |62[ ) and Eq. (|5T1), we find that when translated into the language of JKLW, the results of 
Ref. give Eq. (|5^ ) as the mean fluctuation of the induced chromoelectric fleld. The essential 
nonlinearity of Eq. (H3) is, therefore, absent in this calculation. 



So far, we have only considered the real part of the JKLW kernel, x- Of course, the same method 
can be applied to flnd what is the form of the virtual part a, Eq. (|4l|), that arises from the 
calculation of Ref. |7| . To reproduce the virtual part it is clearly necessary to keep the quadratic 
terms in the relation between and a*. Thus the quadratic terms in Eq. ( |8a| ) will be important in 
this calculation. Other than that, the calculation is straightforward. Again the gauge invariance 
ensures that all the "kinematical" factors of a of Refs. [p|-|l2| are reproduced in the projectile 
gauge calculation and the only difference comes from the difference in the normalization of the 
eigenfunctions. It is clear, therefore, that the result of such a calculation is again the lowest order 
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expansion of a in powers of 7+ 

'W'^i){xi.\^\xi_) 









ab \ 



fbcd{x±\ 

D 



■ I 1 1 1 



D2 



x,y)fbcde''^{y. 



ab 



1 



1 



1 



cd 



1 



1 



cd 



\y±)- (65) 



3.4 The Doubly Logarithmic Limit. 

Before exploring the relationship between the two approaches further in the next section, we want 
to make a comment about the form Eq. (^4|). Although this equation certainly gives x iii general 
as a nonlinear function of the background field a*, this nonlinearity disappears in the double 
logarithmic limit. Following Ref. ||l^ we take the double logarithmic limit as the limit when the 
background field a* does not depend on In this limit, the covariant and the simple derivatives 
commute and it is easy to see that Eq. (|4|) reduces to 



Q2 £)2 



or 



trx = 4tr 



a 



Q2 ■ 



(66) 



(67) 



When substituted into the evolution equation Eq. (|4^) this gives the simple linear double logarith- 
mic DGLAP evolution for the gluon distribution function G oc tra^ (see Ref. for a detailed 
derivation). This is in contrast with the situation discussed in Ref. where the double logarith- 
mic limit of Eq. ( ^3|) was studied. It was shown there that the nonlinearities in Eq. (^) survive in 
the doubly logarithmic limit and, in fact, lead even in this limit to the "almost saturation" of the 
gluon distribution. 

The absence of the nonlinearities is in contradiction with the explicit calculation of Mueller and 



Qiu 1 25] who showed that the QCD evolution of the gluon distribution in the doubly logarithmic 
approximation does indeed contain contributions from higher twist operators. This again under- 
scores our observation that the nonlinearities included in the evolution of Ref. |^ are not the whole 
story. Those are the "kinematical" nonlinearities in the sense discussed in the previous section 
and do not include interesting dynamical effects which come into play when the parton density 
becomes large. 

In fact, the triviality of the doubly logarithmic limit of the calculation of the previous section is 
easy to understand using the intuition of the dipole model approach. In the projectile evolution 
picture, the doubly logarithmic limit is achieved by assuming that, in every step in the evolution, 
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the extra gluon that is emitted into the virtual photon wave function has the smallest transverse 
momentum or, in the coordinate space, has the largest transverse coordinate ||8|. The part of 
the wave function that contains this gluon, therefore, essentially describes one adjoint dipole of 
large transverse size. One leg of this dipole is the newly emitted gluon, while the other leg is the 
remainder partons which are closely bunched together in the coordinate space^. The cross section 
for the scattering of the adjoint dipole in the large Nc limit is simply related to the cross section for 
the fundamental dipole cXadj = 2(Tfund — o"fund ■ Thus, if initially one starts (like in Ref. or Ref. 
from a fundamental dipole, the appearance of the large adjoint dipole in the wave function leads to 
a nonlinear GLR type quadratic term |^] in the evolution equation for the scattering cross section. 
This result indeed has been derived in Ref. |^. However, if we want to consider the evolution of 
the gluon distribution itself, the initial state should contain an adjoint rather than a fundamental 
dipole. This can be achieved by considering "DIS" of a virtual particle that couples to trF^ |27|. 
In this case, in any step in the doubly logarithmic evolution, the state contains only one adjoint 
dipole. The probability for the appearance of a larger dipole in the approximation of independent 
emissions is itself proportional to the number of gluons. The evolution of the gluon distribution in 
this approximation is, therefore, necessarily linear and is merely the simple DGLAP DLA. 



4 More on the target versus projectile gauge 

The discussion of the previous section may seem a little paradoxical on the purely technical level. 
Indeed, we have been calculating the same physical quantity in two different ways. The quantity 
in question is the equal (lightcone) time propagator of the transverse components of the vector 
potential in the target lightcone gauge. The first way of performing the calculation is to work 
entirely in the target gauge as was done in Refs. [^12|. This gives the result Eq. (^). The second 
way to calculate the same quantity is to first calculate the propagator of a* in the projectile gauge 



and then gauge transform the result into the target gauge using Eqs. (^,61). This results in an 
inequivalent expression Eq. (p^). 

4.1 The ie complication. 

Our first aim in this section is to resolve this technical paradox. To do this let us consider in more 
detail the calculation of Refs. p|-p!^ and its transformation into the projectile gauge. The equal 
time propagator of the transverse components of vector potential is calculated in the following 
way (9-12]. One starts with the quadratic part of the action for the small fluctuations of A^. 



Integrating A it is reduced to a quadratic action for the small fluctuations of the transverse 



^These partons are in the adjoint representation of the color group, since together with the extra emitted gluon 
the state must be an overall singlet. 
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components of the vector potential 

S = j (fx(fyA\x)G-^'\x,y)A^{y). (68) 

For the purpose of this discussion we use somewhat simpUfied notations and omit the color indices 
on the fields. The explicit form of G^^ is given in Refs. P-p!^. One then finds properly normalized 
eigenfunctions of 

I d^xA^^p. = '^(A - X')5{p- -p-')5\r - r') . (69) 

Using the complete set of eigenfunctions one constructs the propagator with the standard ie pre- 
scription as 

The limit = y"*", and x~, y~ — > cxd is then taken to calculate x*'^- Clearly, the equal time limit 
selects the on-shell eigenfunctions A = and, therefore, when transforming into the projectile gauge 
it is important to keep track of the ie prescription. The simplest way to do this is to include the 
ie term directly in the action 

St = j (fxd^yA\x) \g-^'\x, y) + ie5'^ 5{x - y)] {y) . (71) 

The propagator Eq. ( [70[ ) is then just the inverse of the quadratic form in Eq. (71) without any 
additional regulators. 

To transform this expression into the projectile gauge one has to use Eqs. (57,31). The gauge 
invariance of the QCD action ensures that the first term in Eq. (^) under this transformation 
transforms into 

a\x)D\x,y)a\y), (72) 

which is exactly the action used in Sec. |2|to calculate the projectile gauge propagator. However, the 

ie term is not so simple. If the transformation Eq. ( |6l|) was unitary, the norm of the field A^ would 

be preserved and the ie term in Eq. (^) would transform into the standard ie f d'^xai{x)ai{x) term 

in the projectile gauge. The problem is that, as we saw in the previous section, the transformation 

Eq. (^Tj) is not unitary. A normalized function a* is transformed into a function A^ normalized not 

to unity but rather to an eigenvalue of the operator O in Eq. (|52|).P| This was precisely the root 

^The fact that the transformation between the two gauges is non-unitary is not unusual. Even though it is a 
gauge transformation and therefore formally unitary, the gauge parameter itself depends on the dynamical field. 
Such transformations are generically non-unitary and do not preserve the scalar product. 
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for the discrepancy between x X- The resulting projectile gauge action can be written asf^ 



<fxaHx)D'^aHx) + 



+ ie I d X 



a\x)a\x) + 2a\f)-)[d' 
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(73) 



Thus, the standard ie prescription in the target gauge is equivalent to a fairly complicated mo- 
mentum dependent prescription in the projectile gauge. Since the calculation of Sec. |2|, following 
Ref. 0, was performed using the standard ie prescription in the projectile gauge the result is, 
indeed, expected to differ from that of Refs. • 



While the technical reason for the difference between the results of Refs. |^-12| and Ref. 0] is clear, 
the physics behind it is not so obvious. In the rest of this section, we will make an attempt to 
understand the physical reason for this difference. 

As we have just explained, the calculation of Refs. |P-12| is equivalent to a calculation in the 
projectile gauge with a nonstandard momentum dependent ie prescription. It is well known that 
such a change of prescription is equivalent to a calculation not in the vacuum state but rather in 
a state which contains gluons [28,^,^]. We, therefore, ask ourselves why the projectile gauge 
calculation should be performed in a state which, on top of the background field ft"*", also contains 
additional gluons. 



4.2 Evolution as Renormalization Group in the projectile gauge. 

To answer this question let us first try to reformulate the projectile gauge calculation of Sec. ^ in 

terms of the Wilson renormalization group akin to the approach of Refs. [pHl^]. The hadron is 

represented as a statistical ensemble of the static fields of the form Eq. (|^) with a statistical 

weight Z[b]. Evolution in x generates induced vector potential which changes the statistical weight. 

Strictly speaking the induced vector potential is not static. It has components in the frequency 

range < A oc However, as long as the frequency of the components of the projectile wave 

function are large enough, one can treat the induced potential as static during the interaction with 

the qq pair. Also, as long as the wavelength of the projectile in the x~ direction is large enough (p+ 

is small) the induced vector potential can be approximated by a S{x~) shaped function. One can, 

equivalently, describe the hadron by a statistical ensemble of V and with some statistical weight 

Z[V,V^. It is clear that both descriptions encode exactly the same information. For this purpose 

one has to define analogs of x ^^id a, that is the (connected) fluctuation correlation functions of 

^"in writing this expression we have made use of the fact that the ie term is important only for functions a' that 
satisfy D^a^ = 0. 
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order as- The resulting evolution equation 
d 



din 



-Z[U, c/t] 



a. 



1 



j—[ZXqq{u,v)\ + 



<5C/t('u)(5C/(f) 



+ 



2\5U{u)6W{v 



is the analog of Eq. (p]). 

Using the formulae of Sec. (|2|) and the appendix we find 



(74) 



(75) 



(x - z)_L • (y - z)^ 



27r2 7 (x — 2:)^(y — z)^ 

(f/(^x)-f/ny±))"'(f/n^±)-^(^±))'^ 
+ (C/(xx) • C/^(^±))"'(^^(y±) • ^(^±))^^ 
- 5'''{u\y^) ■ U{x^)y^ - (C/(xx) • U\y^)y'5^^ 



(76) 



vgq y-^^.u^j - 27r2 7 " {x - z)l{y - z)\ 

'{Uiz^)Y'{Uiy^)-U\z^)-Uix^)y^ 
+ {U{x^)-UHz^)-U{y^)Y\u{z^)y^ 

- {U{x^r\u{y^f' - {U{y^r\u{x^)y'] (77) 

aq is obtained from Ug by replacing C/s by C/^s. The same is true for Xqq and Xgg- Xgg is obtained 
from Xqq by swapping (x, a/5) with (y, 7^) as this exchanges the q and q lines there. We stress 
that Eqs.(|7^^) contain all the information that is contained in the BK equation as well as in the 
equations for higher correlation functions of U that appear in Ref. 

4.3 The snag. 

There is one implicit assumption in this procedure: namely that the a"*" component of the vector 
potential is the only relevant one. One assumes that if, for example, an a* component is generated 



25 



in the evolution it does not affect tlie subsequent evolution of the physical cross section. This is, 
however, not quite right. What is important for the interaction with the projectile is not merely 
but rather the F~^^ component of the color electric field. The interaction between the projectile 
and the target is due to the term in the QCD Lagrangian. The F"* component is the 

Weiszaker- Williams field of the qq pair, while the F"*"* component is generated by the color charges 
in the target. In the eikonal approximation it is true that F~^^ = d^b'^ and, therefore, the coupling 
can be written as b'^J~, where J~ = However, if there is a contribution to F"*"* coming 

from the transverse component of the vector potential, it should be taken into account. 

It is easy to see that such a contribution is indeed generated by the low x evolution. Suppose one 
starts the evolution initially with the background field configuration as in Eq. (|l|). In the first step 
of the evolution one generates both the increment in o"*" and the increment in a*. The two are 
related by the condition Eq. ( |60| ) 

a+ = ^a^ . (78) 

Naively, one would expect that since all the fiuctuation fields in this step have small frequencies, it 
should be true that <C a* and, therefore, it should be safe to forget about a*. The reason this is 
incorrect is that the on-shell solutions for a* are discontinuous at x~ = 0. Therefore, even though 
the field a* is indeed small, it has a large derivative with respect to x~ which contributes to the 
field strength. In fact, the induced chromoelectric field is 

= d'a+ - D+a' . (79) 

Recalling that on-shell a* satisfy the second of the equations Eq. (|60|), we see that the second 
term in this expression is 2B^o} and is of the same order as the first term ^^aK Clearly, even if 
one starts initially from a background which only contains b'^ after long enough evolution a large 
transverse component of the vector potential is generated. When the contribution of the transverse 
component to the field strength is comparable to the contribution of the "+" component, the eikonal 
approximation breaks down and the evolution discussed in Sec. ^ ceases to be valid. It looks indeed 
very natural that in order to take into account the presence of the (potentially large) transverse 
field, the calculation in the projectile gauge should be performed around a state that contains 
transverse gluons apart from the 6"*" background. 

One could try to argue that the transverse part of the vector potential can be somehow gauged 
away and the calculation could still be performed consistently around a pure b'^ background. Even 
if this is possible the evolution of the background defined by such a procedure will be different 
from the evolution of Sec. ^. In any case, we do not see how such "regauging" is possible. 

It is instructive to see in more detail how the gauge fixing works in both the projectile and the 
target gauges and why the two seem to have different status as far as the renormalization group 
structure is concerned. As we mentioned above, the chromoelectric field is created by the color 
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charges in the target. In fact, the whole renormahzation group procedure can be formulated in 
terms of the color charge density rather than the vector potentials themselves, which was in 
fact originally done in Refs. p|-p^. The background vector potentials are found as static solutions 
of classical equations of motion in the presence of the color charge density j"*" = p6{x~) 

F'^ = , 

D'F+' = D'[D%+ - d+b']= j+ . (80) 

An important property of these equations is that for a given p they have infinitely many solutions. 
By considering an arbitrary unitary matrix V{x±,x~) it is straightforward to see that all the 
following are solutions 

b+ = + D'd+b']. (81) 

The difference between the target and the projectile gauges at this point becomes important. In 
the target gauge = 0, the equations reduce to 

B' = iV^d'V, 

D'^d+B' = -i+ . (82) 

We, therefore, get rid of almost all the solutions, the only residual degeneracy being the value 
of the matrix V aX x~ ^ — oo. The imposition of the residual gauge condition d'^A^{x~ — > — oo) 
then removes all solutions except one. In the projectile gauge the situation is very different. The 
condition a~ = does not eliminate any of the infinite number of solutions Eq. (|8l|). The choice 
of the residual gauge fixing is thus crucial to eliminate the redundant solutions. If those are not 
eliminated, the perturbative calculation will be plagued with zero mode problems. The calculation 
in Sec. ^ was in fact performed with the residual gauge fixing d^a^{p~ = 0) = 0. This gauge fixing 
does indeed eliminate all the solutions except the one which has vanishing a* and has therefore 
precisely the form of Eq. (||). 

Now, consider the renormahzation group calculation. Here, we have to integrate out modes which 
have higher frequency p~. In the target gauge this is straightforward: the residual gauge condition 
does not care about frequency. It therefore eliminates nonzero frequency fiuctuation modes which 
do not vanish at — > — oo in the same way as it eliminated the static background solutions with 
this behavior. As a result, the fiuctuation modes have a very similar structure to the background 
field and the induced field is similar to the background. It is, therefore, straightforward to for- 
mulate a selfsimilar renormahzation group transformation in this gauge. The situation is quite 
different in the projectile gauge. The residual gauge condition, although it fixes unambiguously 
the background, has nothing to say about the fiuctuations — it only fixes the static modes! It is 
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impossible therefore to ensure that the fluctuations will have the same form as the static back- 
ground. In fact, as we have seen above, it will not be the case. The first equation of Eq. (^) in the 
projectile gauge is just one of the equations of motion (with or without the external source). This 
means that it is possible to have nonvanishing with vanishing a' only at exactly zero frequency. 
At any finite frequency nonvanishing, a* is required. As we have seen, this a* contributes to the 
induced chromoelectric field, or equivalently to the induced color charge density. The calculation 
thus explicitly lacks a selfsimilar structure and proper renormalization group setup does not seem 
possibl^^ unless extra eikonal approximation is invoked. 

The discussion of this section leads us to conclude that the projectile gauge calculation, as formu- 
lated in Ref. 0] and Sec. § is only valid as long as the eikonal approximation is applicable. When 
the evolution is continued for a large span of 1/x, the eikonal approximation breaks down and 



the higher nonlinear corrections of Refs. P-12| should become important. This is not to say that 



one cannot learn much from this simplified evolution. Quite to the contrary — clearly there is a 
range of x values where this evolution captures the relevant physics. This is particularly true when 
the target is large — the case of a large nucleus discussed in Ref. p. In this case, the eikonal 
cross section is significantly different from the simple perturbative one which assumes single scat- 
tering. The nonlinearity of the evolution becomes important much faster than for a small hadron. 
One, therefore, expects unitarization to appear already within the eikonal regime. Subsequent 
appearance of other nonlinear corrections will not change the fact that the total cross section has 
unitarized. It must, however, affect other more exclusive properties of the process such as the 
structure of final states. The spectrum of the target gauge fields is presumably directly related to 
the spectrum of the emitted gluons |27,3^. Thus, when the evolution of these fields changes even 
locally one expects this change to be visible in the spectrum of final state gluons. Assuming the 
local parton-hadron duality this then has to be mirrored in the spectrum of final state hadrons. 



5 Conclusions 

In recent years several approaches to the evolution of dense gluonic systems in the saturation 
regime were developed. The approaches differ from one another in many technical respects and 
relationship between the physics is also not always clear. In this paper, our aim was to relate two 
of these approaches and thereby to try and reduce the entropy in the field. We have shown that 
the nonlinear JKLW equation of Refs. [^-12] coincides with the BK evolution equation derived 



in Ref. and as long as the gluon field induced by the evolution is small. We have argued 
that the approach of Ref. should break down when the field is large enough so that the eikonal 



^^It may be possible to reformulate RG so that it would include also transverse background fields or equivalently 
finite number of gluons in addition to 6^. This seems, however, to be quite a complicated problem and is far beyond 
the scope of our discussion here. 
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approximation intrinsic in the derivation of Ref. Q ceases to be vahd. We have also argued that 
the evolution of Refs. |^12| when translated into the language of Ref. corresponds to taking 
into account some non-eikonal contributions. 

We should note that our discussion puts into perspective the discrepancy between the double 
logarithmic limit (DLL) of the evolution of Refs. |^12] and the evolution suggested in Refs. [31-33|. 



The AGL equation |^-33| has been shown to arise from the BK equation in the regime where 
the evolution on the projectile side is dominated by production of small size dipoles, or equivalently 
large transverse momentum gluons |8|,|l^. This is a natural regime when the target is a small object 
rather than a large object of the typical hadronic size. It was suggested in Ref. that these 
configurations also dominate in usual DIS in the saturation regime. It seems to us that this point 
warrants further study. In any case this is not the standard DLL, where the evolution on the 
projectile side is dominated by large dipoles. There is, therefore, no reason to expect that the DLL 
of Refs. P-p!^ has much to do with the AGL equation. In fact, as we have shown the DLL of the 
BK evolution is itself extremely simple when considered as the evolution of the gluon distribution 
operator rather than the physical DIS cross section. It turns out to be entirely devoid of nonlinear 
corrections and coincides with the standard DGLAP double logarithmic equation. The DIS cross 
section still evolves nonlinearly and in fact saturates in this limit due to the nonlinear Glauber type 



relation between the cross section and the gluon distribution On the other hand, the DLL 

of Refs. [p|-p!^ is also nonlinear for the gluon distribution and as a result the evolution is slowed 
down already on the level of the gluon distribution |3^] . 

We hope that this paper clarifies to some extent the relationship between the different approaches 
to the nonlinear low x evolution. There are still many questions to be answered. In particular, 
it is desirable to find a more explicit relation between the nonlinearities of the JKLW equation 
and the breakdown of the eikonal approximation and to better understand the physics of these 
nonlinearities. Perhaps the most interesting question concerns the effect of these nonlinearities on 
the structure of the final states. Some work on the analytic understanding of quantities less inclusive 
than the total cross section has appeared recently [p5|-37|. There is also an ongoing numerical effort 



in connection with the heavy ion physics |3S,^| in the framework of the McLerran-Venugopalan 



model [pn,|2T|. Further progress in this direction is extremely important both for our understanding 



of the nonlinear physics and for disentangling linear and nonlinear effects in the existing data. 
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A Appendix 



In this appendix, we give a more detailed derivation of the evolution equation of Sec. |2| including 
the calculation of the contributions where the glue does not interact with the background (target). 

We start with the quadratic action for small fluctuations in the projectile gauge 

1 



S=-{a+[-{d-f]a+-2{d^at){d-a^,) 



ab 



The equation of motion for is 



Substituting this in Eq. ([8^) we get 



h ' 



where 



lab 



ab 



(83) 

(84) 
(85) 
(86) 



is the same as the inverse propagator for a charged scalar field in the presence of a background 
field 6'*'. It has been calculated, for example, in Ref. [16|. The most useful form for us is 

dp' 



1 






-J 

ab J 



2p-(27r)3 



[B{x- - y-)9{p-) - 9{y- - x-)9{-p-)] x 



with p~^ 



2p- 



2p- 



and q =p ■ The color matrix tj^^{x ,y ,z±) is0 



U-,Hx-,y-,z^) = {9{-x-)e{-y-) + e{x-)e{y-))6ab 

+ d{-x-)e{y-)Uab{zi_) + e{x-)d{-y-)ul{z^) 



The on-shell two-point correlator of a+ can be written as 

(a+(x+ = 0,x^,x )a,j"(y+ = 0,y_L,y )> = (^a^(x+ = 0, j;^, j; )4(y+ = 0, y_L, y )^) 



dv 



-dldl 



dp 



1 1 



p- 2(27r) [p-y 



[e{x~ - y-)e{p-) - e{y- - x-)9{-p-)] X 



d'^z± 



(2vr)2" 



(2vr)2 



^■^More rigorously, the structure of is given by U~i^{z±) = e*'*'^ )-*(!/ ))>'(^±) ^ However, the difference between 
this expression and that given in Eq.(p^ only shows up if it is multiplied by derivatives or 5{x~ ) factors. Since 
we encounter no such factors in our calculation we will be using Eq.(|l7]) throughout. 
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We now need to expand the eikonal factors 

V{x±) = Vexp 



dx {b~^ + ga'^){x±,x ) 



(90) 



to second order in the fields. Recahing that the background part of the field b'^ (x 6{x ) and that 
the fiuctuation field a+ is nonsingular at = 0, this becomes 



V{x+ = 0,x±) =Vexp 



/ dx 

J —<x> 



{x'^ = 0,x±,x ) 



X U{x±)Vexp 



r+oo 
Jo ' 



dx a {x =0,x±,x ) 



where 



U{x±) = Vexp 



— I 



/ + 0O 
-oo 



dx = 0, xx, x ) 



is the classical part of the eikonal factor. 
To second order in we have 

V{x& = 0,x_l) = 



<l — ig dx a'^{x~^ = 0,x±,x ) 
J —oo 

— 9^ dx~dy~9{y~ — x~)a~^{x^ = 0,x±,x~)a'^{x^ = 0,x±,y~)> 

J —oo ^ 

X U (xj^) <1 — ig dx a'^{x'^ = 0, x±,x ) 

r+oo . 

— g^ dx~dy~9{y~ — x~)a'^{x'^ = 0, x±, x~)a'^{x'^ = 0, x±, y~) \ , 



or 



V{x^) = U{x^) - 

— ig^ J dw~a'^{x^ = 0,x±,w~)U(x±) + U{x±) J dw~a'^{x~^ = 0,x±,w~)^ 
— g^\ / dw~dz~9{z~ — w~)a'^{x^ = 0,x±,w^)a^{x^ = 0,x±,z^)U{x±) 

J —oo 

r+oo 

+ U{x±) / dw~dz~9{z~ — w~)a'^{x'^ = 0,x±,w'')a'^{x'^ = 0,x±, z~) 
Jo 

+ / dw a'^{x'^ = 0,x±,w )U{x± 



(91) 



(92) 



(93) 



+ -n ^. ..-\TTf^.\ 1^ dz-a+{x+ =0,x±,z-)Y (94) 



— oo 







31 



Similarly, 



V\yL) = U\y,_) + 

+ i9{uHy±)J dz-a+{y+ = 0,y±,z~) + dz' a+ {y+ = 0,y±, z-)U\y±)^ 
-g^\u''{y±) / dw~dz~9{w~ - z~)a'^{y^ = 0,y±,w~)a'^{y^ = 0,y±,z~) 

r+oo 

+ / dw^dz'9{w^ - z')a^{y^ = 0,y±,w^)a^{y^ = 0,y±,z~)U^y±) 







+ 







dw-a+{y+ = 0,y^,w-)UHy±) dz-a+{y+ = 0,y±, z')}. (95) 



Rather than directly calculating the eikonal cross section, we will first calculate the tensor product 
of two eikonal factors and later take the trace over the color indices. We use the notation A® B = 
A^^^B"/^ . To second order in the fluctuation field we have 







/ dw a+(x+ = 0,x_L,'w )U{x^) i^U^y^) / dz a'^{y~^ = 0,y±, z ) 

J —oo J —oo 

rO r+oo 

dw^ a'^ {x^ = 0, x±,w~)U{x±) ® / dz^a'^{y^ = 0,y±,z^)U^{y±) 

'-00 Jo 

/•+00 i-O 

+ U{xi_) I dw^a'^ix^ = 0,x±,w~) (S)U^{y±) / dz^a'^{y^ = 0,y±, z^) 

Jo J -oo 

f + OO f + OO 

+ U{x±) / dw^a'^{x^ = 0,x±,w'') <S / dz^a'^{y^ = 0,y±,z'^)U\y±) 
Jo Jo 



dw dz 9{z — w {x'^ = ^,x±,w )a'''(x^ = 0, z )U {x i_) ® U\y ±) 

r+oo 

dw^a'^{x^ = 0,x±,w^)U{x±) / dz^a'^{x^ = 0, x±, z^) (ii^ U^y±) 







+ 00 



U{x±) / dw dz 6{z —w )a'^{x^ = 0,x±,w )a'^{x^ = 0,x±, z )<S^U^{y±] 











U{x±) ^U\y±) / dw dz 6{w -z )a'^ {y^ = 0,y±,w )a+(y+ = 0, y^, z ) 



oo 

+00 />0 



U{x^)(S) / dw a+{y+ = 0,y^,w )U\y±) / dz a+{y+ = 0,y±,z ) 

Jo J -oo 

r+oo 

U{x±) (g) / dw^dz^e^w' - z")a+(y+ = 0, y±, w")a+(y+ = 0, y±,z-)U\y±)) . (96) 







For calculational purposes it is fruitful to separate the different contributions into the set terms 
where the glue interacts with the background field (the first four terms in Eq. ( |96| ) and the set of 
contributions where there is no such interaction (the other terms in Eq. (|9^)). 
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Let us first consider the contribution where the glue is exchanged between the quark in the negative 
half plane {x~ < 0) in the amplitude and an antiquark in the positive halfplane in the complex 
conjugate amplitude: 

r+oo 

dw- / dz-{at{x+ = 0,x±,w-)a+iy+ = 0,y^,z-)) ■ [(taUix^)) ^ {tbU^y^))] . (97) 

-oo Jo 

One notices that only the 0{—p~) term in Eq. ( |89| ) survives and that, Eq. (^), the color matrix 
Uabiw^,z~,z±) reduces to Uab{z±). 

r+oo 

dw' / dz'{a^{x^ = 0,x±,w~)a^{y+ = 0,y^,z')) = 



oo 



xM dw-e'^v-"^ \ dz-e^'^p-^ j ■Uab{z±). (98) 



The w and z integrations are easily performed 



_,£x:^- , _ +iSJ^z- {2p- 



Noting that 





dw'e "^p-"" I dz-e'^'^p-^ = - ^"'l \ . (99) 
P±^q±^ 



^f.-p^f.-J^feir). (100, 

P J -oo P J P J P 



we have 



TTJ P 

The transverse momenta integrals yield 

72 



Taking the derivatives 



i^'/o — /"'^ (x-.)ifa-.)i 

We use the following identity valid for SU (N) group 

Uab{z^){taU{x^)) {tbU\y±)) = 

= 2tT[taUiz^)hU\z^)] itaUix^)r^it,U\y^)r' 

= ^[N,{Uiz^) ■ Wiy^)y\u\z^) ■ Uix^)y^ - Uix^r^U^ (y^y'] . (104) 
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This contribution, the second term in Eq. 



therefore is 



1 



^TT^ J V 



{x-z)\{y-z)\ 



2N, 



Nc{U{z^)-UHy^)y\uHz^)-U{x^)y' -Uix^rf'uHy^r' ■ (105) 



7/3 



The third contribution, due to the exchange of the gluon between the quark in the positive half 
plane and the antiquark in the negative halfplane, is calculated similarly with the only difference 
that now U~f^{w~ , z~ , zx_) = ?7^^(2;_l) and we pick up the 9{p~) term in the propagator. 



+ 00 



dw 



dz-{at{x+ = Q,x^,w-)a+{y+ = 0,y^,z-)> • [{U{xi_)ta) [U^yM] 



1 

4tt^ 



+ 00 



dp 



P 
1 



d^z^ 



{x - z)_L • {y - z)±_ 
{x-z)\{y-z)\ 



N,{U{x^).U\z^)Y\u\y^)-U{z^)y^ -Uix^r^U^y^y' . (106) 



For the quark self energy correction everything is the same as in Eq. ( 105 ) except that the transverse 
coordinates of the two fields coincide (y^ — > x_l = x_l) 



dw 



dz (a+(x+ = 0,x_L,u; )al [y'^ = ^,x^ = x^, z )) ■ [{taU{xi_)tb) ®U\yA_)\ 



•dp- 



1 

4vr^7 P 
1 



d'z^- 



(x - z)l 

N,tT{{u{x^) ■ t/^(^x))(f/(z±))"^ - {u{x^)y'']uHy±r' . 



(107) 



Finally, the antiquark self energy correction 



+ 00 



dw 



dz {at{x+ = 0,y±,w )a+ {y+ = 0,y^ = y±, z )) ■ [U{x^) {tatbU\y±))] 



1 dp 

47^3 



-oo P 
1 



d^z±- 



2Nr, 



{y-z)l 

iVctr(C/t(y^) • U{z^)) {U^z^))'' - {U^y^))''] U{x^r^ . 



(108) 



Now let us look at the other set of contributions. Since for all of these there is no interaction with 
the background we have, Eq. (|88|) , U'l^^ {w~ , z~ , z±) = 5ab- The z_\_ integral is trivial for all terms 
in this set (it yields a 5 function for the transverse momenta). However, as we want to combine 
both sets of contributions in the final result, this integral will not be performed. 



Consider the first term in Eq. ( pq ) — the exchange between the quark and the antiquark in the 
negative half-plane — 

/ dw-dz-{at{x+ = 0, yx, w')a+{y+ = 0,y^= y±, z')) ■ [(taUix^)) [U^yM] . (109) 
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In the correlator part now both 9 functions survive, corresponding to the two possible orderings of 

w~ and z~ 



dw dz (a+(2;+ = 0,x_L,u; )a^(y+ = 0, y_L, ^ )) 

1 

dp- I 1 



/ p- 2(27r) {p-y 



2 / d p± , ,(x^-z I ) 

^ ' (27r)2'' 



(2vr)2 



I — -^w +1 z 



dw dz [9{w -z )6{p )-6{z -w )9{-p )]e ^■^p-'" e^''2p-" ] ■ d^b 



The w and z integrals now have to be performed with more care. Doing this we get 



-didi 

TT y .1 p 



dp f .2. f d'^P± _+ip^(^^_^^) f d'^q± 



d^z± 



9{p-){- 



1 



(2^)2 

1 



(27r)2 



:)+0{-p-){ 



1 



1 



'■p±'^q_L^ Pi_'^{q±^ - Pi_'^)' P±'^qi.'^ q±'^iqi.^ - P±^) 

Recalling Eq. (p.OOD we obtain for the momentum denominators 



■6. 



ab 



dp 



9{p'){ 



1 



1 



-^)+0{-P'){ 



1 



P~ 
+^dp~ 



1 



p±'^q±'^ 



+ 



1 



-y^^) 



) 



p p±'^q±'^ 



(110) 



(111) 



(112) 



The correlator part is then 



dw dz (a+(x+ = 0,x_L,M^ {y^ = 0,y±, z )) 



+°°dpf 
P' 



{x - z)± ■ {y - z)_L 
{x - z)\{y - z)\ 



■Sab- (113) 



The color algebra raises no problems 



6ab{taU{x±))^{U\y^)tb) = 

= 6ab{taU{x^)r^{U\yM^' 



1 



(114) 



2Nr 



NX'{U\y^) ■ C/(x^))^^ - Uix^r^U^iy^r' 



and so the first term in Eq. ( P6|) is 

1 f^^^dp- 



4vr^ JO P 



{x-z)\{y-z)\ 



2N,V 



NcS'^^U^y^) ■ U{x^)y^ - Uix^r^U^y^r' 



(115) 



35 







The quark to antiquark exchange in the positive half-plane gives 

f +00 

dw-dz-{at{x+ = 0,y±,w-)a-l;{y+ = 0,y±y±,z-)) ■ [iU{x±)ta) ® itbU^iy±))] 

+00 

1 f dp' f ^2^ (x - z)± ■ {y - z)± ^ 
^ {x-z)\{y-z)l 

, aS 



47r^ J p 



2Nr 



Ncy\U{x^).U\y^)Y-U{x^r^U\y^V' . (116) 



Now we combine the two terms that give corrections to the quark line — the fifth and the seventh 
terms in Eq. (p6|). It is easy to see that they have the same color structure and will also yield the 
same transverse structure. 

The color algebra is trivial 

5ah{tahU{x^)) ® U\y^) = 5a,{tahU{x^)r^{U\y^)r' = C/(x^)"%t(y^)75 . (117) 

The fifth term is the same as Eq. (|115| ), but with y_L — > = x_\_ and only the 6{—p~) term 
surviving 

rO 

dw'dz'9{z' — w')(^a'^ {x^ = 0, x±,w')a'^ {x^ = 0, ■ SabitatbU ix±)) U^iy±) = 
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d^z± 



cfp±_ 
(27r)2 



+ 



Ip^ql qliql-pl) 
And the seventh term is the same as Eq. (p,li 



{N^c - 1) 



U{x^)''(^U\yi_V^ . (118) 



dw dz e{z -w )(a+(x+ = 0,a;_L,u; )a'^ {x^ = 0, x±, z ))-6ab{U{x±)tatb)0U\y±) 



vr " 



P~ 



d^z^ 



d'^P± +iv I {x I -z I ) [ d^q± ^+iq,{z, -xi) 
(2vr)2 



p\qi p\{qi-p\) 



(2vr)2 

2Nr 



Uix^)'"^U\y±y^ . (119) 



Adding the two terms Eqs. ( 118| , 119| ) and the performing the transverse integrations we get 



1 



+00 



dp 



47r^ J Q p 



d^z±^-^^-^$^u{x^r''uHy^r' . 



(120) 



■ (x - z)l 2N, 

The correction to the antiquark line — the eighth and the tenth terms in Eq. (96) give similarly 

'dp' >• 1 '^t2 



1 

47^3 



P' 



■{y-z)i 2N, 



(121) 
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Finally, combining all the terms together we get 

, aS , 



{x - z)± ■ (y - z)± 



tr(?7(xx) • U\z±)){U{z^)f^ - N,U{x±r^]uHy±y' 
U{x^r''\tv{uHy±) ■ U{z^)){U^{z^)f - N,U\y^r' 



{x-z)\{y-z)\ 
1 



(x - z)l 

1 



(122) 



This coincides with the result of Ref. Q . When comparing this evolution equation with the results 
of Ref. 0, one should keep in mind that there the evolution is considered with respect to the 



variable C,. The relation between the two evolution equations is given by i (. . . ) 



(...). 



Now taking trace over the color indices we obtain the evolution equation for the scattering cross 
section given in Sec. 0. 
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